Signatures of nonlinear cavity optomechanics in the weak coupling regime 
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We identify signatures of the intrinsic nonlinear interaction between light and mechanical motion in cavity 
optomechanical systems. These signatures are observable in the resolved sideband regime even when the cavity 
linewidth exceeds the optomechanical coupling rate. A strong laser drive red-detuned by twice the mechanical 
frequency from the cavity resonance frequency makes two-phonon processes resonant, which leads to a non- 
linear version of optomechanically induced transparency. This effect provides a new method of measuring the 
average phonon number of the mechanical oscillator. Furthermore, we show that if the strong laser drive is 
detuned by half the mechanical frequency, optomechanically induced transparency also occurs due to resonant 
two-photon processes. The cavity response to a second probe drive is in this case nonlinear in the probe power. 
These effects should be observable with optomechanical coupling strengths that have already been realized in 
experiments. 
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Introduction. Spectacular advances in the quality factor of 
nano- and micro-mechanical oscillators and their rapidly in- 
creasing coupling to optical and microwave resonators have 
given rise to remarkable progress in the field of cavity op- 
tomechanics HIE). This has enabled cooling of mechanical 
oscillators to their motional quantum ground-state [3,4] and 
observations of optomechanically induced transparency [5-71, 
quantum zero-point motion |8] [9), as well as squeezed light 
and radiation pressure shot noise lT0HT2ll . 

The interaction between light and mechanical motion due to 
radiation pressure is intrinsically nonlinear. While several the- 
oretical studies of the single-photon strong-coupling regime 
have been reported recently [ 1 3 - 20 ] , most realizations of cav- 
ity optomechanics are still in the weak coupling limit where 
the coupling rate is much smaller than the cavity linewidth. 
Experiments to date have relied on strong optical driving, 
which enhances the coupling at the expense of making the ef- 
fective interaction linear. Realizations that show promise for 
entering the strong coupling regime include the use of cold 
atoms [10], superconducting circuits [6|, microtoroids [21], or 
silicon-based optomechanical crystals |4|. In the latter, a ratio 
between the coupling rate and the cavity linewidth of 0.007 
has been reported El . and improvements seem feasible. To 
enter the nonlinear regime of cavity optomechanics is of great 
interest, since it is only then that the internal dynamics can 
lead to true non-classical behaviour 1 23 1. 



In this article, we study corrections to linearized optome- 
chanics and identify signatures of the intrinsic nonlinear cou- 
pling that are observable even with a relatively weak optome- 
chanical coupling. The nonlinear effects we discuss come 
about due to the presence of a strong optical drive. We 
show that if this drive is detuned by twice the mechanical 
frequency from the cavity resonance frequency, two-phonon 
processes become resonant. This gives rise to a nonlinear 
version of optomechanically induced transparency (OMIT). 



OMIT has been well studied in linearized optomechanics [24) 
and is analogous to electromagnetic ally induced transparency 
in atomic systems. We point out that the two-phonon induced 
OMIT enables a precise measurement of the effective aver- 
age phonon number of the mechanical oscillator. This pro- 
vides an alternative to sideband thermometry ||8l I9l l25ll26l . 
Furthermore, we show that OMIT also occurs if the drive is 
detuned by half the mechanical frequency due to two-photon 
resonances, and the cavity response to a second probe drive 
is then nonlinear in probe power. We expect these effects to 
be observable for coupling strengths that have already been 
realized in experiments. Their observation would verify the 
intrinsic nonlinearity of the optomechanical interaction and 
thus open up the possibility of generating non-classical states. 

To relate to previous work, we note that a two-phonon in- 
duced transparency |27| can also occur in optomechanical 
systems where the cavity frequency depends quadratically on 
the position of the mechanical oscillator [28]. In addition, the 
effect of ordinary linear OMIT on higher-order optical side- 
bands was studied in Ref. [29 1. 

Model. We consider a standard optomechanical system de- 



scribed by the Hamiltonian H = H sys 
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where a (c) is the photon (phonon) annihilation operator, 
w r (u) m ) the bare cavity (mechanical) resonance frequency, 
and g the single-photon coupling rate. The mechanical po- 
sition operator is X = x zp fZ, where z = c + $ , x zp f = 
-J%/2rmj m is the size of the zero-point fluctuations, and m 
is the effective mass of the oscillator. The constant |a p | 2 is 
included for convenience and simply shifts the equilibrium 
position of the mechanical oscillator. The cavity mode is 
driven by a laser at the frequency cj p . This drive will be re- 
ferred to as the pump, and is described by the Hamiltonian 



H, 



pump 
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FIG. 1. (color online). Depiction of the setup when the pump detun- 
ing is (a) A p ~ — 2a>m and (b) A p ~ — w m /2. Note the difference 
in scale between (a) and (b). A p and A s are the detunings of the 
pump and probe lasers from the bare cavity resonance frequency oj r , 
whereas A p and A s are the detunings from the effective cavity reso- 
nance frequency. 



Let us move to a frame rotating at the pump frequency ui p 
and perform a displacement transformation, such that a(t) — ► 
e -iu) p t ^ _|_ a(t)}. We define A p = u p — w r ^ as the 
pump detuning from the cavity resonance and choose a p — 
iflp/A p . This results in the Hamiltonian H = Hq + H±, 
where 

H Q = -hApa^a + huj m c i c + TiG (c + c 1 ) (a + a f ) , (2) 
Hi = %g (c + c f ) ata. (3) 

We have introduced G = ga p and assumed, without loss of 
generality, that a p is real. The coupling G is enhanced by the 
square root of the average cavity photon number compared 
to g and provides a bilinear coupling between photons and 
phonons. This coupling has been well studied, and it is known 
to give rise to effects such as sideband cooling 01 [30] 1311 
and OMIT Bd El- 

Identifying resonant nonlinear terms. The bilinear Hamil- 
tonian H with A p < simply describes two linearly cou- 
pled harmonic oscillators. By a symplectic transformation, 
we can express Ho in terms of new operators A and C, which 
are annihilation operators for the normal mode excitations of 
the system. These excitations are in general superpositions of 
photonic and phononic degrees of freedom. Up to a constant, 
the Hamiltonian becomes 



H n 



-hA p A^A + hui m c^c . 



(4) 



We will assume that G/oj m <C 1, and that the pump frequency 
cjp does not coincide with the sideband frequencies w r ± w m , 
but rather that |w m ± A p | is on the order of u m . See Fig. fl] 
for an overview of the frequencies involved. In this case, the 
operator A describes excitations that are photon-like, while C 
describes phonon-like excitations. To second order in G/oj m , 
we get A = [l+2A + A_p/(l-p 2 )]a-A + c-A_c t -A + A_pa t 
and C = [l+2A + A_p/(l-p 2 )]c+A + a-A_a t +A + A_ /3 - 1 c+ 
when we define p = oj m /A p and A± = G/(A p ± w m ). The 
normal-mode frequencies are A p = A p (1 — 2A + A_p) and 

Urn = Wm (l + 2A+A_ / 9~ 1 ). 

We can now rewrite the Hamiltonian Hi in terms of the 
normal-mode operators A and C, which results in multiple 
terms. However, since g/n, G/uj m <C 1, we only retain terms 
of nonzero order in G/ujm if they are resonant. We consider 
two different pump detunings. First, if A p ~ — 2a> m , we find 



Hi ties = -%— (A^C 2 + & 2 A 
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This describes processes where one photon-like excitation is 
created and two phonon-like excitations are destroyed, and 
vice versa. On the other hand, if A p ~ — oj m /2, the resonant 
terms are 



Hi* 



, , G 



i t2 (7 + C^A' 
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which describes processes where two photon-like excitations 
are created and one phonon-like excitation is destroyed, and 
vice versa. The Hamiltonian Q combined with ([5J or ([6]) 
gives rise to new effects beyond standard linearized optome- 
chanics, which we study below. 

Equations of motion. We include dissipation and additional 
drives by input-output theory l32l[33ll . The cavity and me- 
chanical energy decay rates are n and 7, respectively. We 
assume that k>7 and that the system is in the resolved side- 
band regime where w m > k. Note that in the presence of 
dissipation, the amplitude a p = £! p x r (A p ), where the cav- 
ity susceptibility is defined as Xr(^) = (k/2 — iuj)~ l . The 
drive strength il p is related to the laser power P p through 
|J7 p | 2 = K cxt Pp/(TiLu p ), where K cxt < n is the decay rate 
of the port through which the cavity couples to the drive. 

The quantum Langevin equations are 11341 



a = — f — — iAp) a - i(G + ga)(c + &) + x^a^ (7) 
c= — ( — 4- iu) m ) c — iG(a + a*) — iga)a + ^7 c in . (8) 

We now introduce a weak second optical drive, the probe, with 
frequency w s close to the cavity resonance frequency w r . We 
define the optical input operator in Eq. |7]) by y / Kd; n (i) = 
e~' l5t n s + ,/K ex t|ext(i) + ^/Kint£int(t), where K int is the in- 
ternal decay rate, K cxt + K; nt = k and 5 — u> s — uj p is the 
frequency difference between the two drives. The frequency 
O s is related to the probe power P s by \£l s \ 2 = K cxt P s /(hiu s ). 
The vacuum noise operators £ oxt obey (£cxt(£)£cxt(^')) = 
S(t - t') and (|L(t)|«t(f )) - and similarly for | int . 
The mechanical oscillator is not driven, but coupled to a 
thermal bath, such that the mechanical input operator obey 
(Cin(t)4(*')) = Kh + l)S(t - *') and (4(t)&,(f)) = 
nth5(t — £'), where n t h is the bath temperature expressed as 
a number of quanta. We will solve Eqs. |7| and <|8j perturba- 
tively in the single-photon coupling g (35). The coupling G 
cannot be treated perturbatively, but we will exploit the fact 
that G/u m < 1. 

The presence of two optical drives gives rise to a beat note 
in the optical intensity at frequency S. This produces an off- 
resonant drive on the mechanical oscillator, as we have as- 
sumed S 7^ ui m . To avoid parametric instability, we must en- 
sure that the cavity frequency modulations due to the coher- 
ent motion induced by this beat note are much smaller than 



the cavity linewidth, i.e. gG\£l s \/(KUi m ) <C n. This is easily 
fulfilled for a weak probe drive with |fi s | < k along with our 
assumptions of G/u m , g/n <C 1. Note that other instabilities 
can also arise [36 1 and must be avoided. 

We can express Eqs. |7]) and ([H} in terms of the normal mode 
operators A and C. This still gives equations with linear cou- 
pling terms whenever dissipation is present. However, in the 
extreme resolved sideband limit K/ui m <C 1, they simplify to 

i = -(|-iA p )i + ^[ J ff l! i] + 7^a in (9) 

d = -(^+iQ m jC+^[H 1 ,C] + ^c in . (10) 



The effective mechanical linewidth is 7 
p 2 



7- 



vn where v 



4A+A_p/(1 — p 2 ) < f or A p < 0. The effective frequencies 



uj m and A p were defined above. Note that \v 



(G/^ m ) 2 « 

1 such that the effective mechanical linewidth is still small 
compared to the cavity linewidth, i.e. 7 <C K. The effective 
mechanical noise operator is defined by y/j c in — ^/7b hn + 
^/n{Xj r ^ + A_£T) when ignoring the beat note and defining 
yfnt; = ^/Kcxt^oxt + ^/Kint^int- Its autocorrelation properties 
are the same as for c; n , but with n t h replaced by the effective 
phonon number n m = (771th + K ^ 2 -)/l- The last term in n m 
describes heating due to radiation pressure shot noise 1 30 , 3 1 1 . 

Two-phonon induced transparency. We start by focusing 
on the case of a pump detuned by twice the mechanical fre- 
quency, A p = — 2o) m , where two-phonon processes are reso- 
nant according to Eq. |5]). Such processes have been studied 
before for systems with so-called quadratic optomechanical 
coupling 11281 . and it has been shown that they can lead to 
OMIT 1 27 1 much in the same way as single-phonon processes 
do with ordinary linear optomechanical coupling [24], We 
will now see that two-phonon induced transparency can also 
occur in the case of linear optomechanical coupling. 

By solving Eqs. |9]l and ( fT0) i perturbatively in the single- 
photon coupling g and transforming back to the original oper- 
ators, we calculate the optical coherence (a(t)} at frequen- 
cies close to the resonance frequency. Defining the probe 
beam detuning by A s = oj s — ui T and the effective detuning 



A s = A s - A p + A p , we find (&(£)) 



-iStz 



where 



a s ,o 1 



2 3 2 Xr (A s )(z 2 ) 



7- i(A s - A p - 2w m ) 



(11) 



and a Sj o = f2 s x r (A s ). We have introduced g\ = gG/oj m , 
which is the effective coupling rate appearing in Eq. Q. The 
first term in Eq. ( fTT| is the response of an empty cavity. The 
second term, given by a = ,g 2 Xr(A s )[(n m + l)Xr(-Wm) + 
n m Xi(ui m )], does not contain a sharp feature and results from 
Raman scattering of probe photons to the sideband frequen- 
cies ui s ± u) m . The last term gives rise to a narrow dip of 
width 27 in the coherent amplitude as well as a group de- 
lay of the input signal. This is analogous to the well-studied 
case of linear OMIT for pump detuning A p = —ui m . In the 



case of A r 



-2(D m , however, the effect is not due to coher- 




ent driving of the mechanical oscillator. The size of the effect 



FIG. 2. (color online). The cavity response ja s | 2 in Eq. {TT} in units 
of no = (2Q b /k) 2 in the case of pump detuning A p = — 2w m , 
The parameters are G/uj m — 0.05, n/cum = 0.1, Q s /k = 0.01, 
n t h — 1, ando; m /7 = 10 5 . Upper solid: g/n — 0.01. Lower solid: 
g/n — 0.03. Dashed: g/n = 0. Dots: Numerical results. 



rather depends on the average mechanical fluctuations through 
(z 2 ) = (2n m + 1). This is also the case for two-phonon in- 
duced transparency with quadratic optomechanical coupling 
ETIl . Note that the average displacement can be increased by 
mechanically driving the oscillator. 

If the system is not in the extreme resolved sideband limit 
K/u m <C 1, additional terms must be included in Eqs. |9]) 
and ( fTO] ). This leads to the same expression as in ( [TT| >, but 
with the changes g\ -> iuj m g 2 G 2 \xr(um)\ 2 Xr{-u m ), and 
(z 2 ) —> (z 2 ) — K/(2iui m ), The latter correction is due to op- 
tomechanical correlations induced by the radiation pressure 
shot noise [37|. In addition, the effective parameters describ- 
ing the mechanical oscillator are also adjusted if K/w m is not 
negligible l38l . 

The cavity response |a s | 2 to the probe drive is plotted in 
Fig.0for g/n = 0.01 and 0.03. The dip in |a s | 2 corresponds 
to a dip in either transmission or reflection of the probe de- 
pending on the experimental setup. 

The result ( fTT| provides a new way of measuring the aver- 
age phonon number of the mechanical oscillator. To see this 
in an easy way, let us assume n/u) m <C 1 and A p = — 2o) m , 
and that the mechanical oscillator is not driven. We de- 
fine the size of the dip d < 1 at A s = as d e 1 — 
/x|a s /a S! o| 2 = 2fiKi (2n m + 1) to lowest order in gi, where 
K\ = 4g 2 /(Kj) is the effective single-photon cooperativity 
and l//x = 1 - 2Rea = 1 - 2 5 2 |xrKi]| 2 (2n m + 1) w 1. 
In the limit where the optical broadening of the mechani 
cal linewidth is significant, i.e. n(G/u>, 
of the dip becomes d = 9(g/K) 2 (2n r[ 
phonon number can easily be measured if g/n is known. This 
could serve as an alternative to regular sideband thermome- 
try [8, 9 25, 26]. Finally, we note that while the two-phonon 
OMIT is a classical effect, its presence in the low-temperature 
limit n m —> is solely due to quantum zero-point fluctuations 
of the mechanical oscillator. 

Two-photon induced transparency. We now consider the 
case of the pump drive detuned by half the mechanical fre- 



2 ^> 7, the size 
- 1) (35) and the 




FIG. 3. (color online). The cavity response |a s | 2 in Eq. {T2} in units 
of no = (20 s /k) 2 in the case of pump detuning A p = — w m /2. 
The parameters are G/oJ m = 0.05, ft/o; m = 0.05, n t h = 0, and 
W m /7 = 10 5 . Upper solid: g/n = 0.1 and Q s /k — 0.4. Lower 
solid: g/n — 0.01 and Q s /k — 15. Dashed: gJK — 0. Dots: 
Numerical results (only available for weak probe drives in our com- 
putations). The slight difference for larger A s results from different 
values of a. 



quency, A p = — d) m /2, giving rise to the Hamiltonian |6|. 
Again, we calculate the optical coherence for frequencies 
close to the cavity resonance frequency, restricting ourselves 
to the regime n/oJ m <C 1. We find (a(t)) = e~ lSt a s with 



a s ,o 1 



2ff 2 2 |Qs,ol 2 Xr(A s ) 
7/2 - 2i(A s - A p - w m /2) 



(12) 



ignoring a very small constant term. We have defined the 
coupling rate appearing in Eq. (|6]l as 32 = 8/3 <?(G/u; m ) 2 . 
There is also an OMIT effect in this case, as seen from 
the last term in Eq. ( [12) , since two probe photons can 
be converted to one phonon and vice versa. We see that 
the cavity response is nonlinear in probe power. The 
dip size for A p = — <D m /2 at A s = becomes d — 
AiiK 2 \2Vl s /k\ 2 = 32[i(g/K) 2 {G/uj m ) 2 \2n s /K\ 2 , where the 
cooperativity is K2 = ^g^li^l) ar, d the second equality as- 
sumes 7^7. The amplitude |a s | 2 for A p = — a) m /2 is 
plotted in Fig. [3] 

Numerics. To corroborate our analytical results, we nu- 
merically solve the quantum master equation (see e.g. Eq. (9) 
in Ref. ifPfl ) using the system Hamiltonian H(t) = Hq + 
H 1 + ih(e~ i5t Ct s a^ — h.c.) with H and H 1 from Eqs. Q and 
(|3). Since the Hamiltonian only contains the single frequency 
5 = uj s — oj p , we can use the continued-fraction method l40l 
to solve for the frequency components of the density ma- 
trix g(t) = X)n=-7V Qne~ mSt , where TV is an integer cut- 
off. From this, we calculate the part of the optical coherence 
(a) rotating at the frequency of the probe drive. Figs. [2] and 
[3] show that the numerical and analytical calculations are in 
good agreement. 

Conclusion. We have studied corrections to linearized op- 
tomechanics and identified signatures of the intrinsic nonlin- 
ear coupling between light and mechanical motion. The sig- 
natures are nonlinear versions of optomechanically induced 
transparency, that come about due to resonant two-photon 



or two-phonon processes in the presence of a strong, off- 
resonant optical drive. These effects are observable even 
when the single-photon coupling rate is smaller than the cavity 
linewidth, and are thus relevant to present day experiments. 

Acknowledgements. We are grateful for financial support 
from The Danish Council for Independent Research under the 
Sapere Aude program (KB), the Swiss National Science Foun- 
dation through the NCCR Quantum Science and Technology 
(AN), the DARPA QuASAR program (JDT), and from the 
NSF under Grant No. DMR-1004406 (SMG). The numerical 
calculations were performed with the Quantum Optics Tool- 
box ATI . 

Note added. During the final stages of this project, we be- 
came aware of related work by Lemonde, Didier, and Clerk. 



[1] M. Aspelmeyer, T J. Kippenberg, and F. Marquardt, 

arXiv: 1303.0733. 
[2] P. Meystre, Ann. der Physik 525, 215 (2013). 
[3] J. Teufel, T. Donner, D. Li, J. W. Harlow, M. S. Allman, K. Ci- 
cak, A. J. Sirois, J. D. Whittaker, K. Lehnert, and R. Simmonds, 
Nature 475, 359(2011). 
[4] J. Chan, T. P. M. Alegre, A. H. Safavi-Naeini, J. T. Hill, 
A. Krause, S. Groblacher, M. Aspelmeyer, and O. Painter, Na- 
ture 478, 89 (2011). 
[5] S. Weis, R. Riviere, S. Delglise, E. Gavartin, O. Arcizet, 
A. Schliesser, and T. J. Kippenberg, Science 330, 1520 (2010). 
[6] J. D. Teufel, D. Li, M. S. Allman, K. Cicak, A. J. Sirois, J. D. 

Whittaker, and R. W. Simmonds, Nature 471, 204 (201 1). 
[7] A. H. Safavi-Naeini, T. P. M. Alegre, J. Chan, M. Eichenfield, 
M. Winger, Q. Lin, J. T. Hill, D. E. Chang, and O. Painter, 
Nature 472, 69(2011). 
[8] A. H. Safavi-Naeini, J. Chan, J. T. Hill, T P. M. Alegre, 
A. Krause, and O. Painter, Phys. Rev. Lett. 108, 033602 (2012). 
[9] N. Brahms, T. Botter, S. Schreppler, D. W. C. Brooks, and 
D. M. Stamper-Kurn, Phys. Rev. Lett. 108, 133601 (2012). 
[10] D. W. C. Brooks, T. Botter, S. Schreppler, T. P. Purdy, 
N. Brahms, and D. M. Stamper-Kurn, Nature 488, 476 (2012). 
[1 1] T P. Purdy, R. W. Peterson, and C. A. Regal, Science 339, 801 

(2013)| 
[12] A. H. Safavi-Naeini, S. Groblacher, J. T Hill, J. Chan, M. As- 
pelmeyer, and O. Painter, arXiv: 1302.6179. 
[13] P. Rabl, Phys. Rev. Lett. 107, 063601 (2011). 
[14] A. Nunnenkamp, K. B0rkje, and S. M. Girvin, Phys. Rev. Lett. 

107,063602(2011). 
[15] A. Nunnenkamp, K. B0rkje, and S. M. Girvin, Phys. Rev. A 

85,051803(2012). 
[16] J. Qian, A. A. Clerk, K. Hammerer, and F. Marquardt, Phys. 

Rev. Lett. 109, 253601 (2012) 
[17] M. Ludwig, A. H. Safavi-Naeini, O. Painter, and F Marquardt, 

Phys . Rev. Lett. 109, 0636 01 (2012)) 
[18] K. Stannigel, P. Komar, S. J. M. Habraken, S. D. Bennett, M. D. 
Lukin, P. Zoller, and P. Rabl, |Phys. Rev. Lett. 109, 013603] 
(2012)| 
[19] A. Kronwald, M. Ludwig, and F. Marquardt, Phys. Rev. A 87, 

013847 (20T3J1 
[20] J.-Q. Liao and C. K. Law. arXiv: 1206.3085. 
[21] E. Verhagen, S. Deleglise, S. Weis, A. Schliesser, and T. J. 
Kippenberg, Nature 482, 63 (2012). 



[22] J. Chan, A. H. Safavi-Naeini, J. T. Hill, S. Meenehan, and 
O. Painter, Appl. Phys. Lett. 101, 081 1 15 (2012). 

[23] Non-classical means steady states where the Wigner function 
has regions of negativity. 

[24] G. S. Agarwal and S. Huang, [Phyl Rev. A 81, 041803 (2010)| 

[25] A. M. Jayich, J. C. Sankey, K. Borkje, D. Lee, C. Yang, M. Un- 
derwood, L. Childress, A. Petrenko, S. M. Girvin, and J. G. E. 
Harris, New J. Phys. 14, 115018 (2012). 

[26] A. H. Safavi-Naeini, J. Chan, J. T. Hill, S. Groblacher, H. Miao, 
Y. Chen, M. Aspelmeyer, and O. Painter, New 1. Phys. 15, 
035007 (2013). 

[27] S. Huang and G S. Agarwal, |Phys. Rev. A 83, 023823 (2011)] 

[28] A. Nunnenkamp, K. B0rkje, J. G E. Harris, and S. M. Girvin, 
Phys. Rev. A 82, 021806 (2010) 

[29] H. Xiong, L.-G. Si, A.-S. Zheng, X. Yang, and Y Wu, Phys. 
Rev. A 86, 013815(2012). 

[30] I. Wilson-Rae, N. Nooshi, W. Zwerger, and T. J. Kippenberg, 
|Phys. Rev. Lett . 99, 093901 (2 007)] 

[31] F. Marquardt , J. P. Chen, A. A . Clerk, and S. M. Girvin, [PhyT| 
Rev. Lett . 99, 093902 (2007)1 

[32] M. J. Collett and C. W. Gardiner, Phys. Rev. A 30, 1386 (1984) 

[33] A. A. Clerk, M. H. Devoret, S. M. Girvin, F. Marquardt, and 
R. J. Schoelkopf, Rev. Mod. Phys. 82, 1155 (2010) 

[34] One should in general take care when including dissipation for 
strongly coupled modes |42|. However, since we consider an 
off-resonant pump and G <C W m , the hybridization is weak and 



the frequencies involved are only slightly renormalized. Thus, 
we need not worry about the frequency shifts induced by the 
strong coupling and can include dissipation in the standard way 
1321 [33 1 without taking the optomechanical coupling into ac- 
count. 

[35] By evaluating the Routh-Hurwitz criterion [43], one finds that 
the unperturbed system is stable when G 2 < ((k/2) 2 + 
Ap)oj m /(4| A p |), assuming A p < and cj m /7 3> 1. The as- 
sumption G/oj m <C 1 ensures that this criterion is satisfied for 
the detunings we discuss here. 

[36] 1. Suh, M. D. Shaw, H. G. LeDuc, A. J. Weinstein, and K. C. 
Schwab, Nano Letters 12, 6260 (2012) 

[37] Technical laser noise will give additional corrections. 

[38] ForA p ~ -2cD m we get v = -8G 2 a^| X rK0| 2 |Xr(3u; m )| 2 , 
u m = w m - 4G 2 a4 1 (3 + ( K /2w m ) 2 )|xr(w m )| 2 |xr(3w m )| 2 , 
and the phonon number is n m = (771th + «:G 2 |xr(3a; m )| 2 )/7. 

[39] Our perturbative approach implicitly assumes that the mechan- 
ically induced cavity frequency fluctuations are much smaller 
than the cavity linewidth, i.e. g^/2n m + 1 <C n. 

[40] H. Risken, The Fokker-Planck Equation (Springer, Berlin, 
1989). 

[41] S. M. Tan, Journal of Optics B: Quantum and Semiclassical Op- 
tics 1, 424 (1999). 

[42] H. J. Carmichael and D. F Walls, Journal of Physics A: Mathe- 
matical, Nuclear and General 6, 1552 (1973). 

[43] E. X. DeJesus and C. Kaufman, Phys. Rev. A 35, 5288 (1987) 



